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Abstract 

The sine-Gordon model is discussed and analyzed within the framework of the renormalization 
group theory. A perturbative renormalization group procedure is carried out through a decomposi- 
tion of the sine-Gordon field in slow and fast modes. An effective slow modes's theory is derived and 
re-scaled to obtain the model's flow equations. The resulting Kosterlitz-Thouless phase diagram is 
obtained and discussed in detail. The theory's gap is estimated in terms of the sine-Gordon model 
paramaters. The mapping between the sine-Gordon model and models for interacting electrons 
in one dimension, such as the g-ology model and Hubbard model, is discussed and the previous 
renormalization group results, obtained for the sine-Gordon model, are thus borrowed to describe 
different aspects of Luttinger liquid systems, such as the nature of its excitations and phase tran- 
sitions. The calculations are carried out in a thorough and pedagogical manner, aiming the reader 
with no previous experience with the sine-Gordon model or the renormalization group approach. 
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I. INTRODUCTION 



The sine-Gordon model was originally proposed as a toy model for inter- 
acting quantum field theories and it has been intensively investigated ever 
since. 

In low dimensional physics, the sine-Gordon model often appears as a 
description of systems with non-quadratic interactions having a strong pining 
effect. Contrary to the quadratic momentum that promotes fluctuations in 
the system, the sine-Gordon potential would like to lock the model field in 
one of the minima of the cosine. The model is particularly useful to describe 
strongly correlated electronic systems in one dimension. 

As it is well know, interacting electrons systems in dimensions higher than 
D = 1 are well described by Landau's Fermi liquid theory. In D = 1, 
however, the Fermi liquid fails due to an instability - known as Peierls in- 
stability - generated by 2kp scattering processes which are particular to one- 
dimensional Fermi "surfaces". In opposition to the Fermi liquid nomencla- 
ture, one-dimensional electronic systems are generically referred to as Lut- 
tinger liquids, after the related work by Luttinger.- 

Throughout the years, many models and formalisms have been proposed 
to describe the special behavior of Luttinger liquids. Of particular interest 
is the bosonization mapping between different Luttinger liquid hamiltonians, 
such as the g-ology model and Hubbard models, and the sine-Gordon model 
for which plenty of anallitic results are available. 
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The S-matrix formalism for the sine-Gordon model and the elementary 
excitations spectrum have been analytically derived- - -. Results for the 
sine-Gordon model form factors^ and finite size correction to the model's 
spectrum^"- are also available in the literature. 

Nevertheless, there still are a number of open questions regarding the sine- 
Gordon model. Exact results for the model's correlation functions are still 
lacking, for example. The evaluation of the model's spectrum has proved 
challenging as well 10 . 

The renormalization group theory is an important analitic tool is this 
context for it provides the understanding of the sine-Gordon model's phase 
transition and the energy scale at which it occurs. 
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II. THE MODEL 



The sine-Gordon model model hamiltonian is given by 



v 



if [II, ip]= j dx -(n' + (d xV Y)-gcos(/3 V ) 



(1) 



where <p = <p(x,t) and II = U(x,t) are canonically conjugated fields, that is: 



d t if(x, t) = vll(x, t) 



(2) 



The model lagrangean writes: 



L[tp] = J dx{dt<p.n - H[n,ip}} 



L[ip] = J dx 
After an integration by parts, the action 



(3) 



S[<p] = J dtL[<p] 



V Q2,„ 1 . o2 



can be written in the following form 
S[(p] = J J dtdx 

where we have assumed that the (/?-field goes to zero at the boundaries of the 
integration plane. 

For a reason that will become clear in Section III.B, it is more convenient 
to work in imaginary time t — > —it, i.e., with the euclidean action iS — > S 
which reads 

S[<p] = S [<p] + Si[<p] 

/I ] 
dx 2<pV 2 x <p with V 2 x = dl + -^d, 

Sj[if} = J dxlj[(f] with lj[ip] = g cos(f3ip) 
where x — > (x, vt) and g = g/v. 



(4) 
(5) 

(6) 



III. RENORMALIZATION GROUP TREATMENT 

A. Conceptual overview on renormalization group theory 

The renormalization group (R.G.) is essentially a theory of scale invariance 
and symmetries. A symmetry or scale operation on a system is a transforma- 
tion that portraits the system's appearance and behavior at different scales 
(where this scale might be a length scale, energy, or any typical scale in 
the system). The system is said to be scale invariant under a certain scale 
transformation when it looks the same at all scales. 

In the R.G. theory, the system's microscopic physical quantities (such as 
mass, charge and interaction parameters) do depend on the scale at which the 
system is observed. The scale invariant properties are the ones that emerge 
from the system's macroscopic structure and are related to the degree of order 
in the system. The global invariant properties are represented by physical 
quantities called order parameters. 

Many states of matter are characterized in terms of scale invariant prop- 
erties and order parameters, e.g.: the crystalline lattice structure in a solid 
that translates into a periodic density of particles, the spins orientation in a 
ferromagnetic material that results in a net spontaneous magnetization, the 
localization of charge in an insulator described in terms of charge density 
waves, and so on. 

When, for some reason, a scale invariant system loses its invariance it is 
said to have undergone a phase transition. The phase transition, i.e. the 
loss of invariance, takes place at a certain critical scale that is typical of each 
system. The critical scale defines an energy, called gap (or mass, in quantum 
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field theory language), that measures the extent of the disturbance in the 
system's order parameter. The phase transition's critical point is defined by 
the values of the system's parameters at the critical scale. 

As an example, picture a perfect solid at T = 0. As the temperature is 
increased, the solid will eventually lose distance invariance at some critical 
length scale that is set by the characteristics of the material. At this scale, 
the lattice correlations cannot compete with the thermal fluctuations and 
the solid structure melts in a fluid. 

Here, we are interested in the scale behavior of the sine-Gordon model (or 
rather of a certain physical system that can be described by the sine-Gordon 
model). 

The next three sections feature a general presentation of the R.G. proce- 
dure. In Sec. B, the decomposition of a generic quantum field theory in slow 
and fast modes is presented; The goal of Sec. C is to express the so-called 
residual action that mixes slow and fast modes in terms of the theory Green's 
function; In Sec. D, a slow modes' effective action is derived through averag- 
ing out the fast modes. The last two section are dedicated to the application 
of the general formalism to the sine-Gordon model. In Sec. E, the model's 
effective action for the slow modes is evaluated. This effective theory is then 
renormalized, resulting in a re-scaled sine-Gordon model. The model's flow 
equation are derived in Sec. F. 
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B. General procedure I - Decomposition in slow and fast modes 



The R.G. procedure as it is presented in this section follows the formu- 
lation by Kenneth Wilson, developed in the late 60's and which awarded 
him the Nobel Prize in 1982. Nowadays, this formulation is routinely called 
"wilsonian approach" to the R.G theory. 

The procedure is based on splitting the theory's field <p(x) in two compo- 
nents corresponding to different momentum-frequency regions of the original 
field's Fourier decomposition. Mathematically, 



where x — > (x, vt), q — > (q,u/v), qx — > (qx + ut), \q\ 2 — > q 2 + u 2 /v 2 and 



with s ~ 1 and where A is a momentum-frequency cutoff. 

Note that, in the original covariant space, the momentum-frequency shell 
would correspond to the unbounded surface between the two hyperbolaes 
q 2 — u 2 /v 2 = and q 2 — u 2 /v 2 = A 2 . Although this surface imposes a cutoff 
in the modulus \q\, individually the coordinates q and u remain boundless. 
Therefore, in the original covariant space, the integration of a function f(q, oS) 
over the shell will diverge if f(q, to) does not decay fast enough. The purpose 
of the imaginary time rotation performed before Eq. (jl]) is exactly to avoid 
complications that might arise from an unbounded shell. 





bulk — \q\ < — , 
s 



shell = — < \q\ < A, 
s 
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In a more compact form, we may write the c^-field as 

<p(x) = v s (x) + 6(fi(x) (7) 

where: 

Jbuik (2i\T 

J shell K Z7r ) 

The (/? s -field contains the so-called slow modes of the original c^-field while 
the 5ip-fie\d contains the fast modes. 

The idea is to obtain the theory's action, written for the c/> field in the full 
momentum-frequency space, in terms of slow and fast mode fields and take its 
average with respect to the fast modes' unperturbed ground state. The result 
of this average is an effective action for the slow modes. A "renormalized" 
theory is thus obtained from the effective one through a scale transformation, 
or renormalization, of the momentum-frequency cutoff. The R.G. statement, 
based on the assumed scale invariance of the theory, is that the original and 
renormalized theories are equal, i.e. that the slow modes' effective theory 
defined in the bulk is equivalent to a scale renormalization of the full original 
theory in the entire momentum-frequency space. This equivalence allows 
the derivation of the theory's R.G. flow equations which comprise the final 
outcome of the R.G. approach. 

Let us proceed by rewriting the action S[ip] in terms of the ip s - and 5ip- 
fields. 
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Since 

f dx5ipV 2 x (p s = [ dx(p s V 2 x 5tp = 



= [**[ j^f «rt q )(-f-£) v We«*«» 

J Jbulk (27r) 2 J shell W \ V 1 ) 



inserting Eq. fl7j) into Eq. (0) gives: 

S M = S [(p s ] + SoM (10) 

The interaction contribution to the action, that is Sj[(f] in Eq. (jBj), will be 
treated via perturbation theory around the slow modes (/Afield. (This pro- 
cedure is similar to the usual saddle-point expansion around a fixed classical 
field configuration.) Up to second order in the perturbation, i.e. to second 
order in the fast modes J^-field, we have 

Si[(p] = Si[(p s ] + J dx a s (x)5(f(x) + J dxdx 5cp(x)b s (x,x / )5ip(x / ) (11) 

where the coefficients a s (x) and b s (x) are given by: 

a - (x) = 6 JM U , (12) 



1 SHjjip] 

2 5(p(x)5ip(x' 



b^y) = L.°lT,J *- (13) 



Substituting Eqs. ( 1TUT) and ( ITTTl into eq. (Hj), the full action S[(p] can be 
written as 

S[<p] = S[<p'\ + 6S[<p a ,6<p] (14) 
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with 

5S[(p s , 5 if] = 5o[^y] + J dx a s (x)5(p(x) + J dxdx 5(p(x)b s (x,x')5(p(x') 

6S[(p 8 ,S(p] = J dxdx' 5cp(x)[5(x-x l ) ^V 2 x ,+b s (x, x')]5(p(x') + J dx a s (x)dcp(x) 

(15) 

and where, in writing the previous equation, we have applied Eq. (JSJ). 

We see from eq. ( TTlj ) that the full action S[(p] splits into two contributions: 
the action S[if s ] for the slow modes and a residual piece 5S[if s ,5if] that 
mixes slow and fast modes. Performing a simple field transformation that 
eliminates the linear term in 5(p(x), the residual action is (at this order in the 
perturbation theory around the slow modes) a quadratic theory for the fast 
modes with a mass- like term given by b s (x, x') that encodes the influence of 
the slow modes as well as that of the interactions. 

In order to derive an effective theory for the slow modes, we will aver- 
age the residual action 5S[ip s ,5ip] with respect to the unperturbed ground 
state of the fast modes 5(p-Re\d operator such that 5S[ip s ,5(p] will become a 
5S e ff[ip s ] and S[ip] will become a S e ff[<p 8 ]- As already pointed out, the R.G. 
procedure corresponds to re-obtaining S[cp] from S e fj[(p s ] via establishing a 
"renormalized" theory Sr[(p] through a scale renormalization of the theory's 
cutoff: A/s -> A. 

The average of the residual action will be more easily evaluated if we 
express 5S[(p s , 5<p] in terms of Green's functions. 
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C. General procedure II - Expressing 5S in terms of Green's functions 

The Green's function Gq(x : x') for the free residual action is defined 
through the equations: 

G- Q \x,x') = 5{x-x')\Vl 

( 16 ) 

/ dx" Gq\x, x")G (x", x') = 5{x - x') 

It follows from the definition that: 

-V 2 x Gq{x, x') = 8{x — x') 

From the above equation we see that Go(x,x') = Go( ') and thus we 
can Fourier transform the equation to write: 



1 T7 2 [ dq r ( n \M x - x ') - [ dq 



Gq(q) = G (q, u) = - 9 < 2 2 (17) 



(2tt) 2 uv ^2 V v 2 J J (2tt) 2 

2 

q 2 + u 2 /V 

The Green's function x') for the full residual action is defined through 
the equations 

G ^(x^ x ) = Gq (x-i — x ) 
^J" dx !! G ^ (x, x") G {x !! ^ x') = ^27r) 2 ^ — "^') 

where E(x,x / ) is the theory's self-energy that accounts for the corrections 
to the free Green's function due to interactions and external fields. 
Substituting the definition (TTfij). it follows that: 



—^ x G(x : x ) — / dx S(x, x ^)G(x , x ) — — r^^^x — x ) 
z J \^) 
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dx 



dq dk dk dq 
(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 



Z(q,k)G(k',q')e 



I l\ iqx+i(k+k')x"+iq'x' _ 



dq 



(&0< 



Jq(x-x') 



dq dq' 
(2tt) 2 (2tt) 2 



dq" 



q")G(-q" ', q^ 



Jqx+iq'x' 



dq dq ' 5(q + q')e iqx+iq ' x ' 



(2tt) 2 (2tt^ 2 



G(g,g , ) = G (^ + g , ) + G'o(9) 



(2^ 



Sf^OGK,!/) (19) 



The previous is the Dyson equation written in terms of the theory's full 
self-energy. In our perturbation theory around the slow modes, developed 
up to second order in the fast modes (which is analogous to second order in 
a saddle-point expansion), the self energy is simply given by: 



E(g, q') = -b s {q, q'), E(x, x') = -b s {x, x') 



(20) 



Thus, Eq. ( 033 ) can be rewritten as: 

G(q : q') = G (q)5(q + q') - G (q) J b s (q : q")G{-q'\ q') (21) 
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Now we can use Eq. ( 12~TT ) to develop an expansion of G(q, q') in powers of 
the interaction coupling constant g. 



* Perturbative expansion of G(q,q') in powers of g 
At zero-th order in g we have: 

G^(q,q f ) = G (q)5(q + q / ) 

Up to first and second orders in g we have, respectively: 



G^(q,q f ) = G (q)5(q + q f )-G (q) 
G^(q,q , ) = G (q)5(q + q f )- 



dq' 



(2tt) 2 
1 



b s (q,q")G (-q")5(-q" + q / ) 



G Q (q)b s {q iq ')G,{-q') 



G^(q,q f ) = Go(q)5(q + q f )-G (q) J ^ «")[ Go(-q")6(-q" + q') 



;G (-q")b s (-q",q')G (-q')} 



G^(q,q') = Go(q)6(q + q l ] 



1 



;G (q) 



(2tt) 2 
dq" 



G (q)b s (q,q')G (-q') + 



(27T) 



b s (q,q")G (-q")b s (-q",q')G (-q') 



And so on. 



Coming back to the residual action 5S[ip s , dip], applying the definitions 
and ({IB|) + Eq. (J2DD) into Eq. (115]), it follows that: 



5S[ip s ,5(p}= / dxdx' S(p(x)G (x, x')5(p(x') + / dx a s (x)5(p(x) 
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Finally, we can achieve a quadratic expression in the fast modes 5ip-field 
through a simple field transformation. Let: 

5cp(x) = (p{x) + r(x) (22) 

Then: 

5S[ip s , 5if] = 5S[(p s , <p,r] = J dxdx' { (p(x)G~ 1 (x, x')(p(x') + 

+ (p{x) [ 2G-\x, x')r{x')+a s {x')6{x-x') 
+ r(x) [ G~ 1 (x, x')r(x')+a s (x')5(x-x') } 

Now if 

J dx' x')r(x') = — -a s (x) 



i.e., 



then: 



r(x) = — - J dx' G(x, x')a s (x) (23) 



6S[(p*,0l = J dxdx'[(p(x)G-\x,x')<p(x') - ^a s (x)G(x,x')a s (x')} (24) 
We are now ready to average 5S[ip s , (p\ in the fast modes' ground state. 

D. General procedure III - Averaging on the fast modes' ground state 



First, let us set up the preliminaries. Note that eqs. (|2"31), ( TTBj) and 
imply that r(x) = r[cp s (x)}. Since, from eqs. (J7j) and 



ip(x) = <-p s {x) + r(x) + </?(x) 

we can redefine the slow modes to incorporate the field r(x) through a trans- 
formation (p s (x) + r(x) —>■ (f s (x) and write: 

(p(x) = (p s {x) + <p(x) (25) 
15 



Now, due to the same argument which led to eq. (TTUj). 



H [U,<p] = J dx V -(H 2 + (d x y) 2 ) 

splits like: 

H \n, if] = H Q [U S , if 3 } + F [fl, 01 (26) 

Let |0)^, |0)^ s and |0}^ be, respectively, the full, slow modes and fast modes 
unperturbed (i?o's) ground states. From eq. (ESj), we have: 

|0)*= |0)^|0)^ (27) 

The goal now is to rewrite the full unperturbed Green's function 

G (x,x') = v(0\<p(x)<p(x')\0)* 

in the fast modes' subspace. So, using Eq. ((23), we can write: 



G (x,x') = ^(0\if s (x)ip s (x')\0y+ *(Q\q>{x)(p(rf)W+ 

+ ^(0\ip s (x)<p(x')\0) lfi + *{0\y(x)y s (x')W 

The last two terms in the previous equation vanish since the field operators 
inside the brackets act on different subspaces. Scale invariance implies that 
the first two terms must be equal, that is, the space-time correlations do not 
depend on the field's momentum-frequency scale. Therefore 

Go(x J x') = 2^(O\0(x)0(x')\Oy 

G o (x,x') = 2^(O\^(O\0(xMx')\Or s \Of 

G (x,x') = 2^(O\0(x)0(x / )\Of (28) 

where, in deriving of the previous equation, we have substituted Eq. ( 1271 ) 
and then used the fact that ^(0|0) ^ = 1. 
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The Fourier transforms of the unperturbed Green's function and of the 
delta-functions are properly redefined in the fast modes' subspace as 

dq 



Go(x, x') = Go(x — x) = / 

J SI 

5(x) = ( 



iqx 



shell (2 71 " 

f dq 

J shell ( 2 7r) 2 " 

S(\q\-A/s) = [ dxe~ iqx 



G (q)e 



iq(x—x') 



(29) 



(30) 



i.e., constrained to the high momentum-frequency shell. 

The effective slow modes' residual contribution to the action, let us call it 
SS e ff[ip s }, is obtained as the fast modes' average of the residual 5S[(p s ,(p\: 



SS eff [^ s }= e<0|fl%',£]|0> 
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(31) 



Substituting Eqs. ( 124} ) and ( 128| ) into Eq. ( 13TT ) we arrive at: 



-Gq(x, x')G 1 {x,x') — -a s (x)G(x,x')a s (x') 



5S e ff[ip s ] = J dxdx' 
Now, from Eqs. ©, ([IS]), (EES), (j2DD and {23), it follows that: 
5S e ff[(p s ] = y dxdx' | -Go(x, x^fGg x') + b s (x,x')\ — -a s (x)G(x, 



/\ S ( l\ 

X )CL [ X j 



5S eff [ip s ] = ^(0|S [<3|0)^+ / dxdx' 



-Gq(x, x')b s (x, x) a s (x)G(x, x')a s {x) 

2 4 

The first term on the right hand side is just a constant and can be absorbed 
through a trivial redefinition of 5S e ff[(p s ], which can be finally written as 



5S e ff[(p s ] = J dxdx' 



-Go(x, x')b s (x, x') a s (x)Go(x, x')a s (x') 

2 4 



(32) 



where, in the last contribution to the integrand, G(x, x') has been replaced 
by Gq(x, x') to keep terms only up to second order in the ^-coupling. 
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If we now replace 5S[(f s ,5(f] in Eq. (COB) by 5S e ff[(p s ] given in Eq. (1321), 
we write down the full slow modes' effective action as: 



S ef f[<p'\ = S[<p a ]+6Seff[<p a ] (33) 

In order to compute the contribution of 5S e ff[ip s ] to the effective theory, 
we need to explicit the dependence of the integrand in Eq. ( 1321 ) on the (p s - 



field which is specific of each particular quantum field theory. The following 
section performs this task for the sine-Gordon model. The ultimate goal is 
to derive the model's re-scaled action. 

E. Application I - The sine-Gordon model re-scaled action 

Let us compute the first term on the right hand side of Eq. (1321): 



Fi[^ s ] = J dxdx'^G Q (x,x')b s (x,x') (34) 

Recalling the redefinition of the theory's Green function in the fast modes' 
subspace Eq. (jZSD, it follows that: 

From Eqs. (EGO) and ©: 

3 2 3 2 
b s (x,x') = -tjgcos([3(p s )5(x - x') = -y/ 7 [(/]£(x - x') 

b s (q jq ')= [ dxdx'b s {x,x')e iqx+iq ' x ' 
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b s (q, q') = Jdx {-^[^ = b s (q + q') (35) 

FiM ^ ta( ^ o) L(l G » (?) 

Applying Eq. (O), 

FxW = -V(q = 0) f dq 1 

J SI 



p2tt p 

F 1 [^] = -b s (q = 0) / / 

J0 J A 



shell (2tt) 2 q 2 + u 2 /v 2 

2ir rA d6d\q\ 1_ 
A As "W M 



F 1 M = -i-6 s (g = 0)ln( S ) 

where in the last step we have used Eq. (1531). 
Finally, using Eq. (Ej), we get 

FiM = -9 (^) y rfar cos(/V) (36) 

where: 

= - ln(s) (37) 

In order to computer the second term in the right hand side of Eq. (1321) , 
we expand a s (x') around to write 

F 2 [if s ] = - J dxdx' ^G (x, x')a s (x)a s (x') (38) 

up to second order in the expansion, as: 
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v 2 f 

= — I dxdx' dtdt G (x — x,t — t) x 

x [(a s {x,t)) 2 + 



+ a s (x, t)d x a s (x, t){x'—x) + a s (x, t)d t a s (x, t){t'-t) + 

+ ^a s (x, t)d 2 x a s {x, t){x'-xf + ^a s (x, t)d 2 a s (x, t){t'-t) 2 + 

+ a s (x,t)d x d t a s (x,t)(x' - x)(t' - t) } 



Or in a simpler form as: 
,2 



F 2 [ip s ] = - — / dxdXdtdTG (X,T) x 

x [(a s (x,t)) 2 



- Xa s (x, t)d x a s (x, t) - Ta s (x, t)d t a s (x, t) + 

+ -X 2 a s (x, t)d 2 x a s {x, t) + -T 2 a s (x, t)d 2 a s (x, t) + 



+ XTa s {x,t)d x d t a s {x,t)} 
Having Eqs. (ESj) and ( 150| ) in mind, we have that: 

G (q,Lu) = v J dXdTG (X,T)e-^ qX+ujT) 



d™d™G (q,uj) = v I dXdT (-i) n+m X n T m G (X,T)e 



-i(qX+coT) 



i(qX+LuT) 



X n T m G (X, T) = - f ^ i n+m d;d™G (q, co)e 

V J shell (27r)^ 

f i n+m f dadu 
v / dXdTX n T m G (X,T) = / -^d^G (q,uj)S(\q\-A/s) = f w {K/s) 

J V J shell K Z7r ) 

(39) 
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Therefore: 



F 2 [<p s ] = ~J dxdt[f d o d o(A/s)(a s (x : t))' 



- fsigo (A/s) a s (x, t)d x a s (x, t) - fgo d i(A/s) a s (x, t)d t a s (x, t) + 
+ (A/s) a s (x, t)d 2 x a s (x, t) + i/flOflg(A/s) a s (x, t)d 2 a s (x, t) + 

+ f d i d i(A/s)a a (x,t)d x dta 8 (x,t)] 

From Eqs. ( TT2I ) and ©, it follows that: 

a s (x) = -^sin(/3^ s ) 

(a s (x)) 2 = (/^) 2 sin 2 (/V) = ^[1 - cos(2/V)] 

a s (x)a^a s (x) = P(J3g) 2 ^(J3<p a )<nB(J3<p')d x , t tp a = ^sm(2^ s )d x ^ s 
a s (x)d 2 ,a s (x) = /3(/^ 

- P V [1 - cos(2/V)] (d,,^) 2 + sin(2/VRV 



2 

a s (x)^a s (x) = -^-[l - cos(2{3tp s )](d x ip s )(d t <p s ) + ^ sin(2^ s )<9,%r' 



Keeping only the non-oscillatory contributions (since the oscillatory ones 
average to zero when integrated in space and time), we have 

(0a) 2 V 



+ — % — f dxdt 



±f d2qd o(A/s)(d x tp s ) 2 + l -f dpi {Als){d t <p s f + f dldi (A/s)(d x <p s )(d t p s ) 



where V = v f J dxdt is the system's volume in space and time. 
Now, from Eq. ®, 



/ W (A/*) = -i / ^d 2 G (q,uj)8(\ q \-A/s) 

v J shell \ Z7T ) 
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/ W (A/») = ~f ^ 2 d*Go(q,u)6(\ q \ - A/») 
%ai(A/s) = — / dqduGoiq, u)5{\q\ - A/s] 

v J shell { Z7T J 



and, from Eq. ( TT7j ) 



v 2 (q 2 + u 2 /v 2 ) 2 \ q 2 + u 2 /v 2 J 

a Ft n ( \ IQquj/v 2 
d q d„Go(q,uj)-- {q2 + uj2/v2)3 

Making a change to polar coordinates: 

2 - v dOd\q\\q\ 4 n A 2fo „ xn . A / . 2s J 



/ W (A/«) = - jT ^ ^gf _L(i _ 4cos 2 (^)) 8(\q\ - A/s) = 



(2tt) 2 |# v v " vm ' ' ttA 3 



/«(A/,) = - jf £ ™ - 4-^)) *( to | - A/.) = 



dOd\q\\q\ 4 . , . 2 ... . \ / \ 1 ^ 



'0 J A/s 

.2tt M 



V 2 7fA 3 



/ a;a; (A/ S ) = / * / Jl C os(fl)sin(e)a(|g| - A/a) = 

jo ja/s l Z7r J ^|y| 

Absorbing the constant term in the previous computation of i 7 ^^] an d 
performing an integration by parts, we have: 



/?Vs 3 /" 

Using Eq. ( 157| ) and expanding around dl = (s = 1), we obtain 

Fs[¥,i= -v8^--8^-jy 
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(40) 



where we have used the fact that the first term in the parenthesis goes to 
zero in the limit of large momentum-frequency cutoff A. 

Substituting the results (j3Sj) and ( HOI ) for the quantities ( 1511 ) and ( 1551 ) into 
eq. ( 1521 gives: 



di3py\ 2 s /^ 2 \ 6 



8ttA 3 



Then, according to Eqs. ( 1551 ) and 
action can be finally written as: 

Seff[<P a ] = dx 



47T 



, the full slow modes' effective 



1 / dl3P 4 g 2 \ ~ s / d//3 2 \ , 



(41) 



To go back to original scale where (p s — > cp we perform the substitutions 

?7 } 



4 



1 



dx — > ~2 / dx; V 2 — > s 2 V 2 



in the above effective slow modes' action. Therefore, the re-scaled action is: 



dx 



(42) 



Under the field transformation (3(p —>■ (p, the re-scaled and original actions 
become: 



dx 



1 + : | <^V 2 ^ + cjs" 2 ( 1 - -f— ) cos(ip) 



2(5 



47rA 3 



47T 



r^V 2 ^+ cjcos(^) 



(43) 
(44) 
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F. Application II - The sine-Gordon model flow equations 

The R.G. statement regarding the theory's scale invariance amounts to 
matching Eqs. ( ]4"5T) and (JMj). As a result of this equivalence one writes 
the model's renormalized parameters in terms of the bare ones and of the 
re-scaling parameter in the following way: 



Defining the differential of a parameter as dX = Xr — X, we can rewrite 
the above equations in differential form as: 



The R.G. flow equations for the sine-Gordon model parameters are given 
in terms of the scale / = — ln(s) + Iq by: 





d(3 2 



3/3V 
47rA 3 



dl 






where: 




(46) 
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IV. KOSTERLITZ-THOULESS PHASE DIAGRAM 



A. Analysis of the flow equations 

A possible solution for the system of first order coupled differential equa- 
tions given by Eqs. ( )4~5"1 ) can be depicted as a path (K(l), u(l)) in the K — u 



plane, where / is a parametric running variable. Each such path describes the 
flow of the initial point (Kq, uq) = (K(Iq), u(Iq)) when the variable / starts at 
/o and runs in the direction of the model's original scale (i.e. to recover the 
full momentum- frequency space). The phase diagram of Eqs. ( 1451 ) is given 



by the collection of all possible paths (K(l),u(l)) in the K — u plane. 

We start the analysis by observing that Eqs. ( ]4~5l ) imply du/dl = dK/dl = 
for u = 0. We say that the phase diagram has a line of fixed points at 
u = meaning that the flow stops when and if it hits that line for some 
/ = /*. In this case, the system parameters will take on the value (K(l*),0) 
for any scale / > /*. In particular, a system with bare parameters (Kq, 0) will 
not flow at all, i.e., a free quadratic model does not get renormalized under 
a scale transformation, as expected. 

Secondly, we see from the first Eq. ( )4~5l ) that the flow of the coupling u: 
(i) points upward inside the K < 1 half of the phase diagram and downward 
inside the K > 1 half if u > 0; (ii) points downward inside the K < 1 half 
of the phase diagram and upward inside the K > 1 half if u < 0; (iii) points 
horizontally at K = 1. 

Therefore, the fixed points along the line u = are unstable for K < 1 
(a point just above or below the segment {K < l,u = 0} will flow away 
from it) and stable for K > 1 (a point just above or below the segment 
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{K > 1, u = 0} will flow towards it). 

The second Eq. ( 1451 ) implies that the flow of the parameter K: (iv) points 
to the left inside the K > half of the phase diagram and (v) points to the 
right inside the K < half, regardless of the value of it, (vi) stops whenever 
K reaches the line K = 0. In this case, according to items (i) and (ii), u 
flows up (vertically) for it > and down (vertically) for it < 0. In particular, 
this shows that no path can possibly cross the K = line where the flow of 
K changes direction. 

Combining the above conclusions, we see that there are three possible 
regions in the phase diagram: 

1) Strong coupling regime: The region of paths (K(l),u(l)) constrained 
to the K < 1 half of the phase diagram and which flow to the regime of 
large \u\ and small \K\. In this regime, the interaction, whose strength is 
proportional to the value of |w|, is said to be relevant. 

2) Vanishing coupling regime: The region of paths (K(l),u(l)) constrained 
to the K > 1 half of the phase diagram and which flow to the regime of 
vanishing u and fixed K. In this regime, the interaction is irrelevant. 

3) Crossover regime: The region of paths (K(l),u(l)) that go from the 
K > 1 into the K < 1 half of the phase diagram, thus initially flowing 
towards a minimum value of \u\ attained at K = 1 (where du/dl = 0). Past 
this point these paths turn into the regime of large \u\ and small positive K. 
In this case, the interaction is said to be marginal. 
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Notice that, according to item (iv) above, there is no region for paths going 
in the opposite direction, i.e. from 0<K<ltoK>l. 

Let us complement this discussion with a simple algebraic analysis. 

We focus on the region around the line K = 1 which is where the interesting 
physics takes place. Then writing 

K = 1 + v (47) 

we can rewrite Eqs. ( 145} ) up to first order in v as: 

W = ~ 2uv 

= -2u 2 v => ^ = -Au 2 v 
= _ U 2 V ^^- = -2u 2 v 

^_ 2l , 2) = J (u2 _ 2(A ._ 1)2) = 

The quantity u 2 — 2(K — l) 2 is an invariant for each solution (K(l),u(l)), 
i.e., 

u 2 {l) - 2{K{1) -l) 2 = c (48) 

where c is a constant (for a given solution) that can be determined, for 
example, by the initial conditions: c = Uq — 2(Kq — l) 2 . 

Now, let (Kf, \u\ — > 0) be the extreme point of a path (K(l),u(l)) that 
flows to or from the line of fixed points u = 0. Eq. ( I4"5j ) implies that: 



a '/ = 1± Vt 
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We see that for Kf to exist we must have c < 0; otherwise the path does 
not flow to or from the line of fixed points. To be more precise, there are 
three possible cases: 



★ c = 

u(l) = ±V2{K{1) - 1) 
fixedpoint — > (Kf,u) = (1,0) 

★ c < 

K0| < V2\K(l) - 1| 



fixedpoints — > (Kf,u) = (1 ± 




★ c > 

\u(l)\>V2\K(l)-l\ 
$ fixedpoints 

Based on the previous qualitative and quatitative analysis, we can draw 
the phase diagram for the sine-Gordon model as in Fig. 1. This is the 
known Kosterlitz-Thouless (K-T) phase diagram. The straight lines u = 
zt\/2(K — 1), given by the condition c = 0, define the boundaries between 
the strong coupling, the vanishing coupling and the crossover regimes. These 
lines are called "separatrices" . The strong and vanishing coupling regimes 
consist of the family of hyperbolas defined by the condition c < (and thus 
"enclosed" by the separatrices) while the crossover regime corresponds to the 
hyperbolas defined by c > ("outside" the separatrices). 
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FIG. 1: Kosterlitz-Thouless phase diagram for the sine-Gordon model. 

B. Gap 

In both the strong coupling and the crossover regimes the flows are towards 
large \u\. At some critical scale in these flows, call it / c , the interaction 
becomes too strong, driving a phase transition in the system. Thus, at / c , the 
system loses scale invariance and the R.G. statement is no longer valid. Based 
on the perturbative nature of the R.G. procedure, a reasonable estimate for 
l c is the scale at which the flow of |it| reaches unity. The system's critical 
correlation length £ c can be assessed through the expression: £ c oc exp(7 c ). 

Since the gap A oc £~ , 

A = exp(- l c ) 

gives an estimate for the gap (except for a multiplicative energy factor) that 
opens up in a system that starts at \uq\ < 1 and flows to the large |w|-regime. 



29 



Our task is to determine the value of / c , and thus of A, as a function of the 
sine-Gordon bare parameters Kq and \uq\ < 1. A first approximation is to 
consider the perturbative R.G. only up to first order in the coupling constant 
g. At this order, we can straightforwardly integrate the flow equation for u 
and determine l c . 

Keeping corrections only up to first order in g would have led to simplified 
flow equations of the form: 

§ = Mi - K) 

f = 

as can be seen directly from Eqs. ( 135T ) and (HSj) by taking u 2 — > 0. 
Since now K = Kq is a constant parameter, we can write: 

c u du' r l 

/ % = 2(1 - K ) / dl' 

u(l) = u exv[2(l-K )(l-l )} 

We see that, for Kq < 1, \u\ increases boundlessly, while for Kq > 1, \u\ 
decreases until it reaches the line of fixed points u = 0. Note that first order 
perturbative R.G. cannot capture crossover paths. In particular, K = 1 rep- 
resents a line of fixed points at this level of approximation. Just to illustrate, 
the sine-Gordon model phase diagram produced by first order perturbative 
R.G. looks as in Fig. 2. 

Coming back to the gap, for Kq < 1 and |uo| < 1 we can write: 

l = \u \exp[2{l- Kq){I c -Iq)} 
/ c = Z + ln(K| 1/2( ^- 1) ) 
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FIG. 2: Phase diagram for the sine-Gordon model produced by first order perturbative R.G. 

Therefore, in first order approximation: 

. A/r ^ v ( coK| 1/2(1 - Ko) K <l(and\u \<l) 
A = A(K ,u Q )={ (49) 

{ AT > 1 

As shown in Fig. 3, for a given \uq\ < 1, the gap decreases with Kq (since 
| wo | < 1) until it reaches zero at the critical value Kq = 1. On the other 
hand, given Kq < 1, the gap increases with \uq\ until it reaches its maximum 
value of Co = exp(— Iq) corresponding to \uq\ = 1. This behavior of the gap 
with Kq and \uq\ is an expression of the fact that the critical correlation 
length decreases as one goes deeper into the strong coupling regime, i.e. as 
| wo | increases and Kq decreases. 

The line K = 1 in Fig. 2 defines the boundary between the gapless (van- 
ishing coupling regime) and gapped (strong coupling regime) regions of the 
phase diagram. The system can undergo a phase transition between the gap- 
less and gapped phases by varying the parameter Kq across the line K = 1. 
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FIG. 3: First order gap as a function of Kq (left) and \uq\ (right). 

A first correction to the first order gap of Eq. ( EES ) and Fig. 3 can be 
achieved by expanding the gapped region into the crossover regime of Fig. 
1 (where marginal paths may start at the region K > 1 but ultimately flow 
into the large \u\ regime). This correction should take into account that the 
boarder line in the K-T phase diagram is no longer given by K = 1, but by 
K = 1 ± u/v2, where the upper sign stands for u > while the lower one is 
for u < 0. 

Based on this and guided by the first-order results of Fig. 3, we can draw 
a qualitative picture for the gap produced by second-order perturbative R.G. 
such as depicted in Fig 4. 

For a given \uq\ < 1, the gap decreases with Kq until it reaches zero at the 
critical value Kq = 1 + litol/v^- The dashed line on the graph represents the 
gap as given by first order R.G. The region between the two curves accounts 
for the contribution of marginal paths to the gap opening. As indicated by 
the second Eq. (145}). for small i^o, the parameter K remains almost constant 
along the second order R.G. flow. In other words, close to the line K = 0, 
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FIG. 4: Qualitative picture for the second-order gap as a function of Kq (left) and |«q| (right). 

the second order flow is essentially vertical like in first order. Therefore, for 
small Kq, the first- and second-order approximations should give roughly the 
same results for the gap, as shown in Fig. 4. As Kq increases, the first and 
second-order gaps depart from each other to die in different critical points. 

Now, given Kq < 1, the dependance of the gap on \uq\ should be similar 
to that of the first order approximation. It is not possible to apply a similar 
qualitative reasoning for the Kq > 1 region because flows starting there 
have a non-trivial marginal behavior. In particular, since the marginal flows 
are "longer", the critical scale l c (at which a gap would open up) might 
exceed the system's cutoff and, in practice, the phase transition might not 
be realizable. The important point anyway is that, deep inside the strong 
coupling sector of the K-T phase diagram where the gap is more relevant 
(larger), A = co|wo| 1//2( ' 1 ~ K °' remains a good quantitative estimate. 
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V. APPLICATIONS IN CONDENSED MATTER PHYSICS 



The main motivation of the sine-Gordon model to condensed matter 
physics is that the model is the bosonized version of the fermionic g-ology or 
Hubbard models for one-dimensional interacting electron systems (the Lut- 
tinger liquids). In this context, the sine-Gordon bare parameters Uq, Kq 
and the non-renormalized velocity v are connected to the original micro- 
scopic couplings defined for the fermionic models. A comprehensive review 
on bosonization methods can be found in the book "Quantum Physics in 
One Dimension", by T. Giamarchi— . 

A. The g-ology model 

The sine-Gordon bare parameters go, Kq and v are related to the ID g- 
ology model's microscopic couplings according to the following expressions 12 






(51) 





where 



Vv - 



(53) 



TTVf 



9v = 



9\\\ ~ 92\\ T 92± 



(54) 



9Av = 9a\\ ± 941. 



(55) 
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and the sub-indexes v = c, s refer, respectively, to the charge and spin sep- 
arated sectors of the full bosonized hamiltonian. In Eqs. ( IMf ) and (151)1). the 
upper signs refer to c and the lower ones to s. 

In the standard g-ology notation, the coupling c/4 corresponds to forward 
scattering between electrons of equal chirality while g^ and g\ correspond, 
respectively, to forward and backscattering between electrons of different chi- 
ralities. Now, the intensity of each such ^-scattering may depend on whether 
the spins of the two interacting electrons are parallel (#||) or anti-parallel 
(<7_i_, in lack of a better notation). Note that for spinless fermions g<i and 
gi processes are identical since one can exchange the outgoing indiscernible 
particles. But once the spin conies in the picture, these two process become 
intrinsically different and contribute to the bosonized theory in different ways, 
as can be seen from the above equations. 

In the general case, when writing models for interacting electrons one is 
concerned with the standard Coulomb repulsion between the particles. In 
the present context, this translates into positive (/-couplings for all processes. 
However, electrons sometimes can interact in an attractive way (as for ex- 
ample, through a phonon mediated coupling). This possibility is taken into 
account by allowing for (some) processes with negative g-couplings. 

From Eqs. (I50f)-(l55|). we see that the Luttinger liquid separates into a 
charge sector described by a model of free bosons with velocity v c and a spin 
sector that maps into a bosonic sine-Gordon model with parameters v s , go s 
and Kq s . Since v c ^ v s charge and spins excitations travel independently in 
the system. 
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From the point of view of the original electronic system, the massless 
charge sector represents electrons in a metallic phase. The behavior of the 
spin sector is not as simple but can be understood in the context of the 
sine-Gordon model phase diagram with bare parameters determined by Eqs. 
(Enj)-(E5D- This phase diagram is depicted in Fig. 5 below. 



-aLl ractPve i nLeracllnns- 




-r«{MU Islve iJiteractiena- 



FIG. 5: Phase diagram for the spin sector of the g-ology model. 

First of all, note that Eq. ( 1511 ) excludes the (K < 0)-half of the full 
K-T phase diagram. In fact, from Eq. (1461). a negative K corresponds to 
an imaginary (5 which, in turn, leads to a hyperbolic cosine in Eq. (CD). 
Although this is certainly a mathematical possibility, it is not the case of 
physical interest. 

Secondly, Eqs. (1qTT)-(154"1) imply that: Kq s > 1 if gu > (for ^-processes 



of comparable intensity), i.e. if go s < (assuming that gn and gi± have 
the same sign), and vice- versa. Therefore, for a repulsive gi-interaction, the 
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only physically meaningful region of the full K-T phase diagram is the one 
bounded above by the line u s = and to the left by the line K s = 1. On the 
other hand, an attractive gi-interaction is described by the region bounded 
below by the line u s = and to the right by the line K s = 1 (and to the left 
by K s = 0). 

Finally, half of the upper crossover regime of the full K-T phase diagram 
was incorporated to the strong coupling regime since, along the remaining 
part of the (now relevant) flows, u s increases monotonically. The lower 
crossover regime, which would have to be "artificially" interrupted at the 
K s = 1 line, can be excluded all over based on the argument of weak inter- 
actions, i.e. small |uo s |. 

The conclusions that can be gleaned from the phase diagram can be sum- 
marized as follows: Repulsive backscattering processes in ID electronic sys- 
tems are irrelevant and the resulting gapless spin excitations behave, in effect, 
as a collection of free bosons that propagate with velocity v s given by Eqs. 
(E2j)-((55D- Attractive backscattering processes flow to the regime of strong in- 
teractions, i.e. are relevant, causing the opening of a gap in the system's spin 
sector. A gapped spin excitation means that the spin (/vfield gets trapped at 
a minima of the cosine and orders, breaking rotational symmetry—. Assum- 
ing that the nature of electronic interactions, i.e. repulsive or attractive, is 
a definite property of a given system, then it is not possible to drive a phase 
transition by varying the pair (Kq s , uq s ) across the point (1, 0). 
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B. The g-ology model at commensurate fillings - umklapp processes 



In ID electron systems with commensurate fillings there is a fourth type 
of interaction known as umpklapp. The correspondent coupling constant is 
termed go, in the g-ology dictionary. The most known is the case of half- 
filling that corresponds to scattering of two left movers to the other side of 
the Fermi level through a momentum transfer of 4kp from the lattice. For 
quarter-filling, an umklapp will be produced by a similar scattering involving 
now four particles with a momentum transfer of 8k p- 

In any case, given that the system is at a commensurate filling, the bare 
parameter go c of Eq. ( 150ft is no longer zero and is associated with a cosine 



perturbation of the type 



+ go c cos(n. y 8hKq c lp — Sx) (56) 

where 



90c = 9 0c = , x 2 , (57) 



[27ray 

n is the order of the commensurability (which affects the amplitude and the 
wave length of the cosine potential) where n = 1 corresponds to half-filling, 
n = 2 to quarter-filling, etc; and the parameter 5 measures the deviation 
(doping) from the commensurate filling. 

The perturbation will oscillate fast due to the phase shift 5x and its space 
integral will vanish unless 5x — > 0. In other words, away from a commen- 
surate filling (finite 5), the umklapp is absent and we recover the previous 
picture of free bosonic charge excitations. But at a commensurate filling 
(5 = 0), the Luttinger liquid separates into two independent sine-Gordon 
models: one for the charge sector with parameters v c , go c and Kq c = n 2 K{) c 

38 



and one for the spin sector with parameters v s , go s and K 0s (with v v , g§ v and 
K 0v given by eqs. (T5QT) - (T5TT) ) . 

Fig. 6 shows the phase diagram for the charge sector of the g-ology model 
at half-filling (n = 1), assuming a positive umklapp coupling g% n . For the 
charge sector, Eqs. (15TT)-(153f) imply that: K 0c > 1 if gn > g 2 \\ + g 2 ±, with 
no implication on the value of go c (that is proportional to #3 ;n and, thus, 
positive). The previous condition leads to a number of possible scenarios. For 
repulsive interactions, it is verified when backscattering (between electrons 
having parallel spins) is more than twice as intense as forward scattering. 
The opposite holds for attractive interactions. If backscattering is repulsive 
and forward scattering is attractive, the condition is always verified. In the 
opposite scenario, the condition is never verified. In general, it is possible to 
drive a phase transition in the charge sector of a ID commensurate electronic 
system by tuning the strength of the interactions so that the pair (Kq c: uq c ) 
moves across the separatrix K c = 1 + u c /V2. In the vanishing coupling 
regime, the umklapp is irrelevant, charge excitations are gapless and the 
system is a metal. In the strong coupling regime, the umklapp becomes 
relevant, the charge excitations develop a gap and the system turns into an 
insulator. In the crossover regime, the umklapp is marginal. 

Another way to drive a metal-insulator phase transition in a ID electronic 
system is by tuning the filling. Given a fixed (Kq c , mq c ) located in the strong 
coupling regime, the system can undergo a metal-insulator phase transition 
by varying 5, i.e. the commensurability parameter. This is a phase transition 
of incommensurate-commensurate type, also known as Mott-transition. 
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FIG. 6: Phase diagram for the charge sector of the g-ology model at half-filling. 



C. The Hubbard model 



The sine-Gordon bare parameters go, Kq and v are related to the Hubbard 
model's microscopic couplings through the equations^ 

v = c 



9o — > 90v — 



-2U 



(27ra) 2 



V = s 



vK V v K 0v = Vf 



(58) 



(59) 



Kq Kq v V 7TOf 



(60) 



where, as before, z/ = c, s and the upper sign refers to c and the lower one to 
s. In the Hubbard model, the coupling U represents an on-site interaction 
of g\ nature with the extra restriction that, since the interaction is local, it 
can only take place between electrons with opposite spins (due to the Pauli 
principle). The Hubbard model can be seem as a simplification on the g-ology 
model. 



40 



From Eqs. (JS3) and ( IfiOj ) and the fact that the i^-parameter is positive, 
Kq s writes in terms of U as: 

K 0s = * (61) 

u_ 



ttvf 

Also here the Luttinger liquid separates into independent charge and spin 
excitations described, respectively, by a free model and a sine-Gordon model 
with their respective parameters. From Eq. (ED), U > Kq s > 1 and vice- 



versa. Therefore, for a system with a repulsive and weak enough Hubbard 
interaction, the pair of bare parameters (Kq s , uo s ) lies inside the u s < 
vanishing coupling, irrelevant, regime of the full K-T phase diagram. In this 
case, only the gapless phase is accessible for the spin system which consists 
of free bosonic excitations. Meanwhile, for an attractive U, the pair of bare 
parameters (Kq s , wq s ) will fall into either the u s > strong coupling regime 
or in the left half of the crossover regime where the interaction becomes 
relevant. In either cases, the spin sector develops a gap and the spin field 
orders. The spin sector's phase diagram is the same as in Fig. 5 obtained in 
the context of for the g-ology model. Here, again, we do not expect a phase 
transition between the gapless and the gapped phases of the spin excitations 
developing in a metal where the nature of the Hubbard on-site interaction is 
either repulsive or attractive. 

If the system is at half-filling, the charge sector develops an umklapp 
interaction of the form 



go. 



cos(\Z87rKo' c ip — 8x) (62) 
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where in the Hubbard model language: 

2U , . 

90c = 77. N2> 63 

For commensurate fillings other than 1/2, the umklapp interaction assumes 
similar expressions. 

From Eqs. ( 15T)1 ) and (I5D1). Kq c is given in terms of U as: 



1 



Oc 



1 + 

TTVF 



(64) 



Thus, at commensurate fillings, the Luttinger liquid separates into two 
independent sine-Gordon models: one for the charge and one for the spin 
sector with their correspondent parameters. Fig. 7 shows the phase diagram 
for the charge sector of the Hubbard model at half-filling. 
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FIG. 7: Phase diagram for the charge sector of the Hubbard model at half-filling. 
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From Eq. (ED, U > =^> K 0c < 1 and vice- versa. Therefore, for a 
repulsive interaction, the pair of bare parameters (Ko c , u$ c ) falls inside either 
the u c > strong coupling regime or in the left half of the crossover regime. 
In both situations the interaction is relevant and the system opens up a 
gap, becoming an insulator. On the other hand, a weak enough attractive 
interaction puts (Kq c , uo c ) inside the u c < vanishing coupling regime where 
the interaction is irrelevant. In this regime, the gapless charge excitations 
remain in the metallic phase. As before, one cannot drive a metal-insulator 
phase transition between the repulsive and attractive portions of the phase 
diagram in a system where the interactions have a definite nature. 

In summary, the Hubbard model describes the following types of ID sys- 
tems of interacting electrons: Away from commensurability, the system is a 
metal described by gapless charge excitations and, if the Hubbard interaction 
is repulsive, gapless spin excitations that preserve rotational symmetry, or 
gapped and symmetry breaking spin excitations if the interaction is attrac- 
tive. For commensurate fillings, the system will be an insulator formed of 
gapped charge excitations and gapless spin excitations for a repulsive inter- 
action, while an attractive interaction leads to a metal with gapless charge 
excitations and gapped spin excitations. 
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